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ABSTRACT 

This  paper  is  concerned  with  the  comparison  of  four  different  non- 
divergent  wind  fields  obtained  from  a  single  geopotential  height  field 
over  a  dense  data  area.      After  developing  the  divergence  equation  of 
the  non-divergent  stream  function,    four  different  stream  functions  are 
obtained  by  modification  and/or  deletion  from  the  basic  equation. 
Isotach  patterns  for  each  strearfi  iuncty>QB  are  computed.     A  comparison 
of  the  four  stream  function  patterns  and  their  corresponding  isotach 
patterns  is  made.      Two  of  the  stream  functions  are  used  to  obtain  the 
non-divergent  wind  field  for  the  initial  wind  conditions  in  a  primitive 
equation  model  developed  by  Krishnamurth      Two  36-hour  forecasts 
are  made.      A  comparison  of  the  forecasts  is  made  at  three-hour 
intervals  up  to  24  hours  and  at  36  hours.     Both  forecasts  are  compared 
to  the  reported  winds  at  12-hour  intervals. 
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11/ .  __.  the  semi- geostrophic  stream  function 
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1.     INTRODUCTION 

It  is  now  evident  that  the  use  of  pure  quasi-geostrophic  equations 
of  motion  in  numerical  forecasting  has  reached  a  plateau.      There  is 
little  evidence  that  they  will  be  much  more  accurate  as  far  as  prognosis 
is  concerned.      W.    Blumen  (1)  has  given  evidence  of  shortcomings  in 
the  quasi-goestrophic  approximations.      These  are  due  to  excessive 
non-linear  "adjustment  times"  before  pure  geostrophic  balance  is 
reached.      This  balance  is  a  requirement  in  all  quasi-geostrophic  models 
but  has  special  significance  in  the  non-linear  baroclinic  models. 
Charney  (2)  implies  that  although  the  quasi-goestrophic  models,   by 
eliminating  the  high  frequency  noise  disturbances,    have  given  us  valuable 
insight  to  the  physical  properties  of  meteorological  significant  waves, 
they  are  only  the  beginning  of  a  new  era  in  numerical  forecasting.     He 
mentions  the  inherent  errors  of  the  geostrophic  approximation  and  that 
further  improvements  in  the  models  would  only  lead  to  highly  complex 
methods  which  would  tend  to  destroy  the  basic  idea  of  the  geostrophic 
approximation  --  its  simplicity.      Charney  suggests  that  the  primitive 
Eularian  equations  as  first  proposed  by  Richardson  are  the  next  step 
in  successful  numerical  forecasting. 

It  is  known  that  extremely  accurate  initial  wind  conditions  are 
needed  in  order  to  successfully  determine  accelerations  and  divergence 
from  the  primitive  equations;  acceleration  and  divergence  are  an  order 


of  magnitude  smaller  than  the  measurable  quantities  in  the  atmosphere 
(3).      Use  of  the  primitive  equations  will  then  depend  on  our  ability  to 
adequately  describe  an  initial  wind  field  that  will:     1)  be  a  close 
representation  of  the  actual  wind,    2)  be  compatible  with  numerical 
calculation  as  far  as  stability  and  time  incrementation  is  concerned, 
and  3)  be  obtainable  from  present  measured  data.      Charney  (2)  develops 
the  "Balance"  equation  which  is  used  to  describe  an  initial  wind  field 
that  meets  the  above  criteria.      (See  Section  2  for  this  development). 

Charney  also  considers  the  geostrophic  wind  components,    u     and 

O 

v     ,    as  part  of  the  initial  wind  conditions  in  the  primitive  equations. 
Through  a  numerical  example  he  shows  that  large  inertiogravitational 
oscillations  occur.     These  results  show  the  existence  and  accumulation 
of  large  amounts  of  divergence,    a  characteristic  of  gravity  waves  but 
not  of  the  meteorological  significant  waves. 

In  this  paper  the  "Balance"  equation  is  developed.     Using  various 
modifications  of  the  equation,    four  different  non-divergent  stream 
functions  are  determined  and  compared  with  their  isotach  patterns. 
Using  Krishnamurti'  s  primitive  equation  model,    two  36-hour  forecasts 
are  made,    one  with  the  initial  winds  being  geostrophic  (a  modification 
of  the  "Balance"  equation),    and  the  other  with  initial  winds  determined 
from  the  unmodified  "Balance"  equation.      The  two  forecasts  will  be 
compared  at  three-hour  intervals  up  to  24  hours  and  at  36  hours.      Both 
forecasts  will  be  compared  to  reported  winds  at  12,    24  and  36  hours. 
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2.     DEVELOPMENT  OF  THE  NON-DIVERGENT  WIND  EQUATIONS 

There  were  two  approaches  to  the  development  of  the  non-divergent 
wind  equation,    henceforth  called  the  "Balance"  equation.      Charney's 
approach  (2)  was  that  of  necessity.      He  found  that  something  better  than 
the  geostrophic  approximation  was  needed  in  order  to  determine  the 
initial  wind  fields  for  use  in  the  primitive  equations.      Examining  the 
actual  divergence  in  the  atmosphere,    he  noted  that  horizontal  divergences 
were  very  small.      Therefore,    he  said,    the  horizontal  wind  components 
can  be  described  accurately  as  partial  derivitives  of  a  non-divergent 
stream  function,  \jj  ,    as  follows: 


"# 


v  = 


t. 


By  placing  these  values  into  the  divergence  equation,    then  eliminating 
insignificant  terms,    Charney  arrives  at  the  "Balance"  equation.     He 
implies  that  since  this  equation  is  a  representation  of  non-divergent 
winds,    and  since  the  inertiogravitational  winds  are  characterized  by 
relatively  large  divergences,    then  the  equation,    if  satisfied,    will 
prohibit  the  occurrance  of  gravity  waves.     He  mentions  collaborators 
who  have  shown  the  same  results  and  proved  his  implications. 

A  more  sophisticated  approach  to  the  development  is  that  of 
Thompson  (5).     His  intention  was  to  re-examine  the  filtering  problem 
of  the  geostrophic  approximation  and  to  see  if  it  could  be  improved  by 
providing  the  filtering  (elimination  of  gravity  waves)  with  less 
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restrictions.      To  do  this  it  was  necessary  to  determine  the  feature  of 
the  equation  of  motion  that  generates  the  gravity  waves,    and  then  to 
eliminate  this  feature  with  some  kind  of  filtering  approximation. 
Thompson  accomplishes  this  by  obtaining  the  frequency  equation  of  the 
linearized  (by  perturbation  method)  divergence  equation.     He  determines 
that  the  roots  of  the  frequency  equation  which  correspond  to  gravity 
waves,    are  those  due  to  the  total  derivative  of  the  divergence  in  the 
divergence  equation.      By  removing  this  total  derivative  there  is 
eliminated  the  possibility  of  having  any  gravity  waves  when  the  modified 
divergence  equation  is   satisfied.     This  modified  equation  is  the  "Balance" 
equation  as  determined  by  Charney. 

From  Helmholtz,   velocity  can  be  expressed  as  the  sum  of  non- 
divergent  and  divergent  parts: 

V  =  kxV^/    +  VX    •  (2-D 

Starting  with  the  horizontal  equation  of  motion 


dV 
dtiH 


f  kxV 


V$      ■  (2-2) 


the  horizontal  divergence  equation  is  obtained: 

a  v-y  +  ufty.v)  +   vg(v-v)  +7u.%v  +  w,av=_7.fkxV.y^. 

at  ox  3  y  o  x  c?y 

(2-3) 
Removing  the  total  derivative  of  horizontal  divergence  from  (2-3)  leaves, 

|v.Vu+J|v.Vv  =  _V.  fK>cV-V2$    •  (2-4) 

ox  dy  z 

Substituting  the  following  expressions  into  equation  (2-4): 

V^kxyi//  u  =  -l^  v=-|J^.        (2-5) 


d  y 
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the  "Balance"  equation  of  the  non-divergent  wind  is  obtained: 


VfV^/  =  v2f>  -  2J 


(2-6) 


This  equation  defines  the  non-divergent  wind  from  a  height  field,    z, 
and  eliminates  inertic-gravitational  waves  as  long  as  it  is  satisfied. 

Using  equation  (2-6)  Krishnamurti  (4)  develops  expressions  for 
the  following  four  stream  functions: 

(1)  l//„,    the  geostrophic  stream  function  by  dropping  the  non-linear 
term  on  the  right  side  of  (2-6)  and  treating  the  coriolis  parameter,    f, 

as  a  constant,    giving 

V2y/  =  iV2f   »'      g-*h- 

(2)  \1/t  ,   the  linear  balanced  stream  function  by  dropping  the  non- 
linear term  and  treating  f  as  a  variable,    giving 

V2y/  =  -f  (V^*Vf)+rV2$  ■ 

(3)  l/Ayi  >  the  semi-geostrophic  balanced  stream  function  by 
evaluating  the  non-linear  term  using  geostrophic  wind  components 
giving 

v2y/ -{(vy/-vo+lv2$  -¥[«,.  *J. 

(4)  l/Z-NTOj   the  complete  balanced  stream  function  by  making  no 
changes  or  assumptions  in  the  "Balance"  equation. 

The  numerical  considerations  of  convergence  in  solving  for  the  com- 
plete balanced  stream  function  are  discussed  by  Krishnamurti. (2-6). 
Generally,    the  requirement  for  a  convergent  solution  is  that  the  equation 
be  elliptic.      This  requirement  is  met  if: 


v  y  >  |  • 
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3.     PROCEDURES 

Using  200-mb  data  for  1200Z,    13  April  1964,    an  accurate 
geopotential  height  field  (50  meter  interval)  was  plotted  within  the 
following  area:    25N  to  60N,    70W  to  135W.     These  contours  represent 
the  geostrophic  stream  function,  U/„.     Krishnamurti  (4)  uses  this 
height  field  to  compute  the  stream  functions,  1//^,   U^Nl'    anc*  ^N2* 
Figure   1  is  the  geopotential  height  field.     Figures  2,    3  and  4  show  the 
stream  functions,  W^,  U/ jaj  and  \J/t^2'    resPectively,   plotted  with  the 
height  contours  on  a  polar  stereographic  projection.     Stream  function 
interval  is  50  X  10     m     sec"    . 

Taking  the  expression  for  the  non-divergent  wind: 

V  -kxV^/  , 

wind  speeds  were  computed  at  over  200  points  in  the  fields  of  \Jf    , 
llf-KTi  and  I/Ant?'    an<^  isotachs  were  drawn.     Isotachs  of  the  geostrophic 
wind  were  drawn  to  speeds  computed  from: 

V  =  k  X  g^z    . 
8  fo 

Figures  5  through  8  are  the  isotach  fields  of  the  four  stream  functions, 

Three  graphs  of  wind  speed  versus  latitude  were  drawn.     The 

first  is  along  the  line  through  the  isotach  maximums,    labled  line  "A"; 

the  second  along  the  95th  meridian,    labeled  line  "B";  the  third  along 

the  115th  meridian,    labeled  line  "C".     These  graphs  are  Figures  10, 

11  and  12  respectively. 
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Reported  winds  are  shown  in  Figure  9.     Figure  13  is  a  table 
comparing  the  reported  winds  to  each  of  the  four  computed  winds. 

4.     DISCUSSION  OF  WIND  FIELD  COMPARISONS 

Examining  the  streamline  patterns,    Figures  2,    3  and  4,    we  see 
that  definite  cross-contour  flow  exists.     Generally  the  flow  is  to  lower 
heights  upstream  of  close  contour  spacing  and  to  higher  heights  where 
contours  spread.      The  linear  balanced  stream  function,    being  closest 
to  the  geostrophic  wind  in  computational  simplicity,    exhibits  the  least 
cross -contour  flow.     The  complete  balanced  stream  function  shows 
the  most  cross-contour  flow  and  therefore  corresponds  more  closely 
to  statistical  observations  of  this  phenomenon  (6).     The  flow  exhibited 
by  the  non-linear  stream  functions  is  a  direct  manifestation  of  the 
theories  on  geostrophic  adjustment  (7). 

Only  the  complete  balanced  stream  function  shows  strong  anti- 
cyclonic  flow  over  New  Orleans.     This  is  the  area  of  the  sub-tropical 
jet  stream.     The  anti-cyclonic  flow  seems  to  be  excessive  but  may  be 
due  to  boundary  values  in  computations  or  a  lack  of  data. 

It  must  be  remembered  that  the  streamlines  represent  the  non- 
divergent  wind.      The  divergent  parts  of  the  wind  are  small.     Krishnamurti 
says  that  calculations  for  this  map  time  give  divergent  winds  that  are 
at  no  place  greater  than  five  percent  of  the  total  wind. 

Streamline  spacing  corresponds  well  to  contour  spacing  (geostrophic 
streamlines).     This  is  reflected  in  the  isotach  patterns  of  Figures  5 
through  8,    and  the  three  graphs,    Figures  10,    11  and  12.     The  linear 
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balanced  winds  are  close  to  geostrophic  winds.     Krishnamurti  has 
shown  (4)  that  at  low  levels  where  strong  cyclonic  and  closed  circula- 
tions exist,   the  linear  balanced  winds  are  far  too  large.     At  200-mb 
the  geostrophic  wind  maximum  in  the  trough  will  approach  the  complete 
balanced  wind  when  corrected  for  curvature.     The  semi-geostrophic 
balanced  winds  are  close  to  the  complete  balanced  winds.     However, 
only  the  complete  balanced  winds  give  a  good  representation  of  the  sub- 
tropcial  jet. 

Figure  13  shows  a  numerical  comparison  of  the  computed  and 
reported  winds.      The  overall  mean  algebraic  difference  indicates  that 
the  balanced  winds  are  closer  to  the  reported  winds  than  the  geostrophic 
winds  are.     A  comparison  of  this  type  at  the  200-mb  level  must  be 
viewed  with  caution.     It  is  known  that  high  values  of  vertical  wind  shear 
exist  at  this  level.     Small  errors  in  balloon  height  can  give  unrepresenta- 
tive reported  winds.     Looking  at  the  isotach  patterns  of  the  computed 
winds, we  also  see  large  horizontal  shears.    In  areas  of  high  winds  and 
large  shears,    the  downwind  position  of  the  balloon  from  the  reporting 
station  could  be  in  a  sector  of  the  computed  wind  isotach  field  that 
contains  a  velocity  range  of  up  to  20  knots.     It  would  be  very  difficult 
to  include  the  turning  of  the  wind  with  height  to  get  an  accurate  balloon 
position  in  order  to  compare  the  reported  and  computed  wind  at  that 
point.     It  might  be  advantageous  to  make  a  smooth  isotach  analysis  of 
the  reported  winds  and  use  this  to  compare  the  computed  winds.     There 
still  will  exist,    however,    the  problem  of  balloon  location  and  which 

winds  have  errors. 
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Lacey  (9)  compared  the  geostrophic  winds  and  the  non<- divergent 
complete  balanced  winds  with  the  reported  winds  at  850-mb,    700-mb, 
500-mb  and  300-mb  over  the  United  States.     His  comparison  was  made 
for  both  direction  and  speed.     He  found  that  on  the  average  the  com- 
plete balanced  winds  were  much  closer  to  the  reported  winds  in  both 
direction  and  speed  when  the  complete  balanced  wind  and  the  geostrophic 
wind  differed  by  more  than  five  degrees  or  five  knots.     He  also 
determined  that  adding  the  divergent  winds  to  the  non-divergent  complete 
balanced  winds  produced  little  change  in  the  non-divergent  wind  pattern. 
The  divergent  winds  were  about  two  to  three  percent  of  the  total  wind. 

From  the  above  discussion  it  can  be  concluded  that  the  complete 
balanced  wind  is  the  most  accurate  description  of  the  true  wind.     The 
question,    "are  the  complete  balanced  winds  accurate  enough  for  use 
as  initial  wind  conditions  in  the  primitive  equations?  "  will  be  discussed 
in  the  next  section. 

5.     THE  NUMERICAL  FORECASTS 

In  the  Introduction  it  was  mentioned  that  Charney  proposed 
balanced  winds  as  sufficient  initial  conditions  for  numerical  prognosis 
using  the  primitive  equations.     His  main  point  was  that  the  geostrophic 
wind  did  not  filter  the  "noise"  while  the  balanced  winds  did.     He  used  a 
simple  model  to  prove  this  point.     However,    Phillips  (8)  has  shown  that 
because  Charney' s  model  was  too  simple,    he  overlooked  the  fact  that 
using  only  the  non-divergent  winds  as  initial  conditions  (  y  *   V  =  0  at 
t  =  0),   the  forecast  using  an  actual  baroclinic  atmospheric  model 
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would  still  contain  considerable  "noise".      Phillips  claims  that  the  noise 
generated  by  the  restriction  to  non-divergence  could  be  eliminated  by- 
including  the  divergence  obtained  from  geostrophic  approximations  in 
the  initial  data. 

Krishnamurti  has  developed  a  primitive  equation  model  in  which 
he  has  incorporated  the  ideas  of  Charney,    Phillips,    Thompson  and 
others.     His  initial  wind  field  is  a  combination  of  the  non-divergent 
wind  determined  from  the  "Balance"  equation  and  the  divergent  wind 
obtained  by  computing  divergence.     Any  combination  of  computing 
methods  to  calculate  balanced  winds  and  vertical  motions  (thus, 
divergence)  can  be  used  to  determine  the  initial  conditions.     He  has 
done  extensive  work  on  partitioning  the  atmosphere  by  using  his  model 
with  initial  conditions  determined  as  follows: 

a.  The  non-divergent  stream  function,  l//-vr?>    by  using  the 
complete  non-linear  balance  equation. 

b.  Vertical  velocity  expressed  in  a  complicated  omega  equation 
involving   Il/Vt?'    an<^  ^he  divergent  velocity  potential,  \. 

c.  The  divergent  velocity  potential,  X  >    expressed  in  the  mani- 
festation of  the  continuity  equation: 

From  the  initial  conditions  obtained  in  the  above  manner,    a  36-hour 
forecast  was  obtained. 

A  much  simpler  procedure  would  be  to  use  the  quasi-geostrophic 
approximation  in  all  phases  of  calculating  the  initial  wind  conditions. 

Though  this  approximation  itself  has  proven  to  be  inadequate  if  used 
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solely  as  the  initial  wind  condition,    one  might  ask  that  if  it  were  com- 
bined with  quasi-geostrophic  divergence  as  suggested  by  Phillips, 
would  it  still  be  inadequate?     A  second  forecast  was  made  using 
Krishnamurti's  model  with  the  initial  conditions  determined  as  follows: 

a.  The  non-divergent  stream  function,  \lf  ,    determined  from  the 

rg 

geostrophic  approximation. 

b.  Vertical  velocity  expressed  in  the  quasi-geostrophic  omega 
equation. 

c.  The  continuity  equation  expressed  as: 


V2*g=-5T- 


w 


In  both  forecasts,    grid  distance  was: 

y  =  2-1/2  degrees  latitude, 

x  =  y  sine  latitude. 
The  area  covered  was  from  55W  to  135W  and  from  25N  to  60N.      Time 
steps  were  ten  minutes.     North  and  south  boundary  conditions  were 
fixed.     The  eastern  boundary  was  extended  eastward  six  grid  distances 
and  combined  cyclically  with  a  western  boundary  extended  six  grid 
distances  westward.     Data  used  were  the  height  field  from  0000Z, 
12  April  1964.     Figures   14,    1'5  and  16  show  the  contours,    geostrophic 
wind  isotachs,    and  reported  wind  isotachs,    respectively. 

6.     DISCUSSION  OF  RESULTS  OF  THE  FORECASTS 

Figures   17  through  25  show  the  complete  balanced  and  quasi-geo- 
strophic component  u  and  v  fields  at  three-hour  intervals  up  to  24  hours  and 
at  36  hours,    respectively.     These  figures  yield  the  following  observations: 
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a.  At  three  hours  (Figure  17)  there  are  already  considerable 
differences  between  the  balanced  and  geostrophic  components. 

In  the  u  field  at  55N,    75W,    there  exists  a  negative  5  m-sec" 
center  in  the  geostrophic  field,    but  a  positive  15  m-sec        in  the  balanced 
field.      The  elongated  center  on  both  over  the  northwestern  United  States 
if  15  m-sec        larger  in  the  balanced  field.      The  50  m-sec"1  center  off 
Washington  in  the  balanced  field  is  only  20  m-sec        in  the  geostrophic 
field  and  is  displaced  15  to  20  degrees  west.     The  sub-tropical  jet  over 
Florida  is  expressed  well  in  the  balanced  field  but  is  displaced  south 
in  the  geostrophic  field.     Many  irregularities  appear  in  the  geostrophic 
field  while  the  balanced  field  exhibits  a  smooth,    realistic  pattern. 

In  the  v  fields  there  is  much  more  agreement.     This  should  be 
expected  because  of  the  fixed  boundary  values.      We  note  that  the 
positive  area  centered  at  105W  on  both  fields  is  much  smaller  in  area 
and  strength  in  the  balanced  field.     The  opposite  is  true  for  the  negative 
area  at  125W.     A  look  at  Figures   15  and  16  reveal  that  the  balanced 
field  is  more  realistic. 

b.  At  six  hours  (Figure  18)  the  differences  noted  at  three  hours 
between  balanced  and  geostrophic  components  are  even  more  pronounced. 
There  are  now  strong  indications  of  inertio- gravitational  oscillations 

in  the  geostrophic  field.      Lookingat  the  u  fields  we  notice  that  a  large 
high  latitude  negative  center  extends  the  length  of  the  chart  in  the  geo- 
strophic field,    while  the  balanced  field  shows  only  a  decrease  in  the 
positive  zonal  flow  in  the  same  area.      The  similarity  over  northwestern 
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United  States  which  existed  at  three  hours  has  almost  disappeared  even 
though  the  balanced  field  shows  little  change  during  the  three  hours. 
The  balanced  field  of  the  v  component  shows  little  change  in  three  hours, 
but  the  geostrophic  v  field  shows  a  great  increase  of  negative  components 
at  higher  latitudes.     The  overall  similarity  that  existed  at  three  hours 
is  disappearing  at  six  hours. 

c.  At  nine  hours  (Figure   19)  the  balanced  components  are 
behaving  well.      There  is  a  slight  easterly  progression  of  the  u  maximum 
at  130W,    45N,    and  the  ridge  line  in  the  v  field  along  the  west  coast  has 
easterly  movement.     In  the  u  field  we  see  the  sub-tropical  jet  maximum 
moving  eastward. 

The  geostrophic  u  field  has  become  strongly  negative  above 
50N.     The  indication  of  a  sub-tropical  jet  is  becoming  distorted,    and 
it  is  becoming  difficult  to  distinguish  the  developing  trough.     The  v 
field  now  shows  good  continuity  with  the  past  three  hours  and  is  becom- 
ing  similar  to  the  balanced  v  field. 

d.  The  balanced  components  show  good  continuity  at  12  hours 
(Figure  20).     The  geostrophic  u  field  is  now  rapidly  losing  the  high 
latitude  negative  center  that  existed  at  nine  hours.     The  geostrophic  v 
field  retains  continuity  but  there  is  evidence  of  oscillations.      The 
negative  center  at  125W  is  lagging  the  balanced  negative  center  by 
about  ten  degrees. 

e.  The  balanced  components  at  15  hours  (Figure  21)  still  have 
good  continuity  and  retain  pattern  shape.      The  geostrophic  v  field  seems 
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to  have  settled  down  with  a  pattern  development  similar  to  the  balanced 
field.      The  geostrophic  u  field  has  lost  its  negative  values  but  still 
shows  oscillations.      This  field  has  only  half  the  positive  velocity  as 
the  balanced  u  field,    and  no  marked  centers  exist. 

f.  For   18  hours  (Figure  22)  and  21  hours  (Figure  23)  the  same 
trends  exist.     The  geostrophic  u  field  has  again  established  a  negative 
area  above  50N  while  the  rest  of  the  field  bears  little  resemblance  to 
the  balanced  one. 

g.  At  24  hours  (Figure  24)  the  geostrophic  u  field  has  the  negative 
area  in  high  latitudes  with  no  indication  of  a  jet  stream  anywhere. 

The  pattern  slightly  resembles  the  balanced  pattern  but  the  values  are 
all  much  lower.      The  balanced  u  field  appears  to  be  developing  strong 
positive  values  in  the  lower  left  corner.     This  may  be  due  to  the  cyclic 
feed  through  from  the  lower  right  maximum  which  existed  at  the  start 
of  the  forecast.      The  v  fields  show  pattern  similarity  but  the  geostrophic 
field  has  a  great  predominance  of  positive  values;  the  negative  areas 
having  small  values. 

h.     The  36-hour  prognostic  u  and  v  fields  (Figure  25)  show  a 
pattern  breakdown.      The  balanced  u  field  has  strong  distortion  from 
cyclic  feed  through,   but  maintains  a  resemblance  of  the  synoptic 
situation  in  its  pattern  shape.     The  geostrophic  u  field  shows  strong 
distortion  with  a  negative  area  undercutting  a  positive  one.     The  balanced 
v  field  has  distortion  in  the  western  half  with  weak  centers  indicating  a 
very  weak  ridge.      The  geostrophic  v  field  shows  a  counter-clockwise 
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rotation  of  the  prevailing  pattern  in  the  western  half,    with  a  very- 
strong  negative  center  just  east  of  Florida. 

It  can  be  concluded  from  the  above  discussion  that  the  quasi- 
geostrophic  approximation,    even  by  including  divergence,    is  unstable 
and  gives  rise  to  large  inertiogravitational  oscillations  when  used  for 
initial  conditions  in  the  primitive  equations.     From  the  limitations  and 
boundary  values  applied  in  this  model,    the  conclusion  is  most  evident 
in  the  u  component  field. 

A  further  comparison  of  the  prognoses  was  made  by  comparing 
the  reported  wind  speeds  to  the  prognostic  winds.      Prognostic  winds 
were  obtained  by  vector  addition  of  the  u  and  v  fields.     Figures  26,    27 
and  28  show  the  reported  winds  and  the  prognostic  balanced  and  geo- 
strophic  winds  at  12,    24  and  36  hours,    respectively.     The  table  of 
Figure  29  shows  definitely  that  the  balanced  winds  are  much  better  than 
the  geostrophic  winds  for  the  12  and  24-hour  prognoses.     At  36  hours 
too  much  distortion  and  cyclic  feed  through  exist  to  make  a  reliable 
comparison,    for  the  geostrophic  prognosis  appears  to  be  better. 

The  major  conclusion  in  this  numercial  forecast  comparison  is 
that  the  balanced  winds  give  the  prognosis  a  wind  pattern  that  resembles 
the  synoptic  situation.     The  geostrophic  winds  gave  rise  to  considerable 
large-scale  "noise"  and  unrealistic  patterns. 
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